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We present measurements of the process pp — > WZ + X — > £ 1 vgiit at y/s — 1.96 TeV, where I and 
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£' are electrons or muons. Using 1 fb 1 of data from the DO experiment, we observe 13 candidates 
with an expected background of 4.5 ± 0.6 events and measure a cross section o(WZ) = 2.7t\'l P b - 
From the number of observed events and the Z boson transverse momentum distribution, we limit 
the trilinear WW Z gauge couplings to -0.17 < Xz < 0.21 (Akz = 0) at the 95% C.L. for a form 
factor scale A = 2 TeV. Further, assuming that Agf = Ak z , we find -0.12 < An z < 0.29 (X z = 0) 
at the 95% C.L. These are the most restrictive limits on the WW Z couplings available to date. 

PACS numbers: 14.70.Fm, 13.40.Em, 13.87.Ce, 14.70. Hp 



The SU(2)l x U(1)y structure of the standard model 
(SM) Lagrangian requires that the massive electroweak 
gauge bosons, the W and Z bosons, interact with one 
another at trilinear and quadrilinear vertices. In the SM, 
the production cross section for pp — ► WZ + X, a(WZ), 
depends on the strength of the WWZ coupling, gwwz — 
— e cot 9\y, where e is the positron charge and 9w is the 
weak mixing angle. At y/s = 1.96 TeV, the SM predicts 
3.68 ± 0.25 pb 0. Any significant deviation from 
this prediction would be evidence for new physics. 

The WWZ interaction can be parameterized by a gen- 
eralized effective Lagrangian 0, HI with CP-conserving 
trilinear gauge coupling parameters (TGCs) gf , Kz, and 
Xz that describe the coupling strengths of the vector 
bosons to the weak field. The TGCs are commonly 
presented as deviations from their SM values, i.e. as 
Agf = gf — 1) Akz = Kz — 1, and Az, where Xz = in 
the SM. Since tree-level unitarity restricts the anomalous 
couplings to their SM values at asymptotically high en- 
ergies, each of the couplings must be parameterized as a 
form factor, e.g. Xz(s) = Xz /(1 + s/A 2 ) 2 ], where A is the 
form factor scale and s is the square of the invariant mass 
of the WZ system. New physics will result in anomalous 
TGCs and an enhancement in the production cross sec- 
tion as well as modifications to the shapes of kinematic 
distributions, such as the W and Z bosons transverse 
momenta. Because the Fermilab Tevatron is the only 
particle accelerator that can produce the charged state 
WZ + A, this measurement provides a unique opportu- 
nity to study the WWZ TGCs without any assumption 
on the values of the WWj couplings. Measurements of 
TGCs using the WW final state 0,11, H, 0, @] are sensi- 
tive to both the WWj and WWZ couplings at the same 
time and must make some assumption as to how they are 
related to each other. 

WZ production measurements and studies of WWZ 
couplings have been presented previously. The DO Col- 
laboration measured awz = 4.5jJ" @ pb, with a 95% C.L. 
upper limit of 13.3 pb, using 0.3 fb -1 of pp collisions at 
y/s = 1.96 TeV Q. The observed number of candidates 
was used to derive the most restrictive available limits 
on anomalous WWZ couplings. More recently, the CDF 
Collaboration measured awz = 5-Ojli | pb using 1.1 fb -1 
of pp collisions at ^fs = 1.96 TeV but did not present 
any results on WWZ couplings. 

This communication describes a significant improve- 
ment to the previous DO analysis. Not only is the data 



sample more than three times larger, but an improved 
technique is used to constrain the WW Z couplings. In- 
stead of merely the total number of observed events, the 
number and the pr distribution of the Z bosons (pf ) pro- 
duced in the collisions are compared to the expectations 
of non-SM WWZ couplings, significantly increasing the 
power of the WWZ coupling measurement over previous 
measurements 0,0]. 

We search for WZ candidate events in final states with 
three charged leptons, referred to as trileptons, produced 
when Z — > £ + £~ and W — ► i'v, where t and t! are e 
or // . SM backgrounds can be suppressed by requiring 
three isolated high-p^ leptons and large missing trans- 
verse energy ($j>) from the neutrino. The combined 
branching fraction for these four possible final states 
(eee, ee/x, [i[ie and /x/x/x) is 1.5% 

DO is a multipurpose detector [13] composed of several 
subdetectors and a fast triggering system. At the cen- 
ter of the detector is a central tracking system, consist- 
ing of a silicon microstrip tracker (SMT) and a central 
fiber tracker (CFT), both located within a 2 T super- 
conducting solenoidal magnet. These detectors are opti- 
mized for tracking and vertexing at pseudorapidities [131 ] 
|ry| < 3 and \r}\ < 2.5, respectively. The liquid-argon and 
uranium calorimeter has a central section (CC) covering 
|ry| < 1.1, and two end calorimeters (EC) that extend 
coverage to \n\ w 4.2, with all three housed in separate 
cryostats [14j. An outer muon system, covering |ry| < 2, 
consists of a layer of tracking detectors and scintillation 
trigger counters in front of 1.8 T iron toroids, followed 
by two similar layers after the toroids (lij . 

Electrons are identified by their distinctive pattern of 
energy deposits in the calorimeter and by the presence 
of a track in the central tracker that can be extrapo- 
lated from the interaction vertex to a cluster of energy 
in the calorimeter. Electrons measured in the CC (EC) 
must have < 1.1 (1.5 < |ry| < 2.5). Electrons must 
have transverse energy Et > 15 GeV and be isolated 
from other energy clusters. A likelihood variable, formed 
from the quality of the electron track and its spatial 
and momentum match to the calorimeter cluster and the 
calorimeter cluster information, is used to discriminate 
electron candidates from instrumental backgrounds. 

Muons tracks are reconstructed using information from 
the muon drift chambers and scintillation detectors and 
must have a matching central track with px > 15 GeV/c. 
Candidate muons are required to be isolated in the 
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FIG. 1: $ T versus dilepton invariant mass of WZ candi- 
date events. The open boxes represent the expected W Z 
signal. The grey boxes represent the sum of the estimated 
backgrounds. The black stars are the data that survive all 
selection criteria. The open circles are data that fail either 
the dilepton invariant mass criterion or have $ T < 20 GeV. 
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FIG. 2: The reconstructed Z boson pr of the WZ candidate 
events used in the WWZ coupling parameter limit setting 
procedure. The solid histogram is the expected sum of signal 
and background for the case of the WWZ coupling param- 
eters set to their SM values. The dotted and double dotted 
histograms are the expected sums of signal and background 
for two different cases of anomalous WWZ coupling parame- 
ter values. The black dots are the data. The final bin is the 
overflow bin. 



calorimeter and tracker to minimize the contribution of 
muons originating from jets [l6j]. 

Events collected from 2002-2006 using single muon, 
single electron, di-electron, and jet triggers were used for 
signal and background studies. The integrated luminosi- 
ties [17] for the eee, eep,, /U/tte, and final states are 
1070 pb" 1 , 1020 pb" 1 , 944 pb^ 1 , and 944 pb" 1 , respec- 
tively. There is a common 6.1% systematic uncertainty 
on the integrated luminosities. 

The WZ event selection requires three reconstructed, 
well-isolated leptons with pt > 15 GeV/c. All three lep- 
tons must be associated with isolated tracks that origi- 
nate from the same collision point and must satisfy the 
electron or muon identification criteria outlined above. 
To select Z bosons, and further suppress background, 
the invariant mass of a like-flavor lepton pair must fall 
within the range 71 to 111 GeV/c 2 for Z — > ee events, 
and 50 to 130 GeV/c 2 for Z — > n/j, events, with the mass 
ranges set by the mass resolution. For eee and /i/x/i decay 
channels, the lepton pair with invariant mass closest to 
that of the Z boson mass are chosen to define the Z boson 
daughter particles. The $ T is required to be greater than 
20 GeV, consistent with the decay of a W boson. The 
transverse recoil of the WZ system, calculated using the 
vector sum of the transverse momenta of the charged lep- 
tons and missing transverse energy, is required to be less 
than 50 GeV/c. This selection reduces the background 
contribution from ti production to a negligible level. 

WZ event detection efficiencies are determined for 



each of the four final states. Monte Carlo (MC) events are 
generated using pythia [18| and a GEANT [l9| detector 
simulation and then processed using the same reconstruc- 
tion chain as the data. Lepton identification efficiences 
are determined from study of Z bosons in the DO data. 
The average efficiencies for detecting an electron or muon 
with E T (p T ) > 15 GeV are (91 ± 2)% and (90 ± 2)%, 
respectively. The trigger efficiency for events with two 
(or more) electrons is estimated to be (99 ±1)%. For 
events with two or three muons, the trigger efficiencies 
are estimated to be (91 ± 5)% and (98 ± 2)%, respec- 
tively. The kinematic and geometric acceptances range 
from 29% for the eee decay mode to 45% for the /i/x/i 
decay mode. It is also necessary to account for r — > e, [i 
final states of W Z that contribute to the signal. The 
number of r events expected to satisfy the selection cri- 
teria is 0.67 ± 0.11 events. These are treated as signal in 
the cross section analysis, but are treated as background 
in the TGC analysis. Table U summarizes the efficiency 
determinations. 

A total of 13 WZ candidate events is found. Figure [1] 
shows ]prp versus the dilepton invariant mass for the back- 
ground, the expected W Z signal, and the data, including 
the candidates. Table U also details the number of candi- 
dates in each channel. 

The main background for WZ -> t'vll are Z + X 
events where X is a jet that has been misidentified as an 
electron or muon. We assess the background from Z+jets 
production by using an inclusive jet data sample that is 
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selected with an independent jet trigger. Events charac- 
teristic of QCD two-jet production are used to measure 
the probability, as a function of jet Et and 77, that a sin- 
gle jet will be misidentified as a muon or electron. Next, 
sub-samples of ee+jets, e/i+jets, and /z/z+jets events are 
selected using the same criteria as for the WZ signal 
except that the requirements for a third lepton in the 
event are dropped. The single jet-lepton misidentifica- 
tion probabilities are then convoluted with the measured 
jet distributions in the dilepton+jets sub-samples to pro- 
vide an estimate of the background from iT+jets events. 
The contribution for all four decay modes totals 1.3±0.1 
events. 

All other backgrounds are determined using MC. Non- 
negligible backgrounds include SM ZZ production, Zj 
production, and W*Z, WZ*, or Wj* production. We 
define these processes as three-lepton final states pro- 
duced through the decay of one on-mass-shell and one 
off-shell vector boson. These backgrounds and their de- 
termination are described as follows. 

ZZ production becomes a background when both Z 
bosons decay to charged leptons and one of the final 
state leptons escapes detection, thus mimicking a neu- 
trino. The total contribution from ZZ production is 
0.70 ± 0.08 events. 

Z7 final states can be misidentified as WZ events if 
the photon is mis-reconstructed as an electron and there 
is sufficient ~$ T . We estimate the 11 + 7 contribution 
using Z ± 7 MC [2(| combined with the probability for 
a photon to be misidentified as an electron (4.2 ± 1.5%) 
determined from studies of events with photons. This 
process is a background only to the eee and /z/ze final 
states. The total contribution is 1.4 ± 0.5 events. 

The contribution to the background from off-shell 
bosons should be nearly the same as occurs in similar 
processes and a fraction relative to the expected signal is 
determined from ZZ MC events generated using PYTHIA. 
It depends on the decay channel and varies from 8% for 
the ee/z mode to 15% for the /z/z/z mode. The uncertain- 
ties include all of those used for the signal plus an addi- 
tional 16% systematic component to account for uncer- 
tainties in the off-shell component of the MC. The total 
contribution of this background is 0.99 ± 0.19 events. 

To cross check the background estimates, we compare 
the number of observed events with that expected when 
we do not apply the dilepton invariant mass selection and 
the $ T selection. We expected to observe 12.5±1.4 events 
from signal and 62.9 ± 8.4 events from backgrounds. We 
observe the 78 events shown in Fig. Q] 

The SM predicts that 9.2±1.0 WZ events are expected 
to be observed in the data sample. The probability for 
the background, 4.5 ± 0.6 events, to fluctuate to 13 or 
more events is 1.2 x 10~ 3 , which translates to a one- 
sided Gaussian significance of 3. Oct, determined by using 
a Poisson distribution for the number of observed events 
in each channel convoluted with a Gaussian to model 



TABLE I: The numbers of candidate events, expected sig- 
nal events, and estimated background events, and the overall 
detection efficiency for the four final states. 



Final 


Number of 


Expected 


Estimated 


Overall 


State 


Candidate 


Signal 


Background 


Efficiency 




Events 


Events 


Events 




eee 


2 


2.3 ±0.2 


1.2 ±0.1 


0.16 ±0.02 


ee/z 


1 


2.2 ±0.2 


0.46 ±0.03 


0.17 ±0.02 




8 


2.2 ±0.3 


2.0 ±0.4 


0.17 ±0.03 


fififi 


2 


2.5 ±0.4 


0.86 ±0.06 


0.21 ±0.03 


Total 


13 


9.2 ± 1.0 


4.5 ±0.6 





TABLE II: One-dimensional 95% C.L. intervals on X z , Agf , 
and Akz for two sets of form factor scale, A. 



A = 1.5 TeV 



A = 2.0 TeV 



-0.18 <\ z < 0.22 
-0.15 < Agf < 0.35 



-0.17 < X z < 0.21 
-0.14 < Agf < 0.34 



-0.14 < Ak z = Agf < 0.31 -0.12 < Ak z = Agf < 0.29 



the systematic uncertainty on the background. A likeli- 
hood method 21[ taking into account correlations among 
systematic uncertainties is used to determine the most 
probable WZ cross section. The cross section a(WZ) 
is 2.7^J'3 Pk> where- the ±lcr uncertainties are the 68% 
C.L. limits from the minimum of the negative log likeli- 
hood. The uncertainty is dominated by the statistics of 
the number of observed events. 

By comparing the measured cross section and pf distri- 
bution to models with anomalous TGCs, we set one- and 
two-dimensional limits on the three CP-conserving cou- 
pling parameters. A comparison of the observed Z boson 
Pt distribution with MC predictions is shown in Fig. [2] 
We use the Hagiwara-Woodside-Zeppenfeld (HWZ) [22j 
leading-order event generator processed with a fast detec- 
tor and event reconstruction simulation to produce events 
with anomalous WW Z couplings and simulate their ef- 
ficiencies and acceptances. The HWZ event generator 
does not account for r final states, and as a result, we 
treat the 0.7 event r contribution as background for the 
WW Z coupling limit setting procedure. The method 
used to determine the coupling limits is described in 
Ref. [23)]. Limits are set on the coupling parameters \z, 
Agf , and Akz- Two-dimensional grids are constructed 
in which the parameters Xz and Agf are allowed to vary 
simultaneously. Table |TT] presents the one-dimensional 
95% C.L. limits on Xz, A<?f and Akz- Figure [3] presents 
the two-dimensional 95% C.L. limits under the assump- 
tion Agf = An z for A = 2 TeV. The form factor 
scale, A (24j, associated with each grid, is chosen such 
that the limits are within the unitarity bound. 

In summary, we present the results of a search for 
WZ production in 1.0 fb^ 1 of pp collisions at y/s = 1.96 
TeV. We observe 13 trilepton candidate events with an 
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FIG. 3: Two-dimensional 95% C.L. contour limit in Agf = 
Akz versus AXz space (inner contour). The form factor scale 
for this contour is A = 2 TeV. The physically allowed region 
(unitarity limit) is bounded by the outer contour. The cross 
hairs are the 95% C.L. one-dimensional limits. 



expected 9.2 ± 1.0 signal events and 4.5 ± 0.6 events 
from background. This gives an observed significance 
of 3. Oct. We measure the WZ production cross section 
to be 2.7^^3 pb, in agreement with the SM prediction. 
We use the measured cross section and distribution to 
improve constraints on WWZ trilinear gauge couplings 
by a factor of two over the previous best results. 
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